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Abstract 

The model-based approach, used by the JeNo code to predict jet noise spectral directivity, is 
described. A linearized form of Lilley’s equation governs the non-causal Green’s function of interest, 
with the non-linear terms on the right hand side identified as the source. A Reynolds-averaged Navier- 
Stokes (RANS) solution yields the required mean flow for the solution of the propagation Green’s 
function in a locally parallel flow. The RANS solution also produces time- and length-scales needed to 
model the non-compact source, the turbulent velocity correlation tensor, with exponential temporal and 
spatial functions. It is shown that while an exact non-causal Green’s function accurately predicts the 
observed shift in the location of the spectrum peak with angle as well as the angularity of sound at low to 
moderate Mach numbers, the polar directivity of radiated sound is not entirely captured by this Green’s 
function at high subsonic and supersonic acoustic Mach numbers. Results presented for unheated jets in 
the Mach number range of 0.51 to 1.8 suggest that near the peak radiation angle of high-speed jets, a 
different source/Green’s function convolution integral may be required in order to capture the peak 
observed directivity of jet noise. A sample Mach 0.90 heated jet is also discussed that highlights the 
requirements for a comprehensive jet noise prediction model. 


1. Nomenclature 

Xj Cartesian coordinate components; i = 1 denotes axial direction 

Vi Cartesian coordinate components 

$ displacement vector components 

£ displacement vector magnitude 

t time 

p density 

p pressure 

T temperature 

T r temperature ratio, plenum to ambient 

Vj velocity component 

«; fluctuating velocity component 

c speed of sound, = yf/RT 

U axial mean velocity 

Uj velocity at jet exit 

Dj jet diameter 

D Green’s function directivity function. 
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L Pridmore-Brown operator 

n reduced pressure \\ = {\ly) In (p/pj 

8jj Kronecker delta 

8 delta function 

y specific heat ratio 

n' dependent variable for pressure 

Q source component, (eq. 3) 

/ frequency 

f m m th mode of solution to second-order compressible Rayleigh operator 
h enthalpy 

r Cartesian coordinate, radial component 

R spherical coordinate component, radius; distance to observer 

6 spherical coordinate component, polar angle 

cp spherical coordinate component, azimuthal angle 

co frequency, 2uf 

M Mach number 

St Strouhal number, 0)DJ{2kUj) 

He Helmholtz number, coD j /( 2 nc,) 

k acoustic wavenumber 

r time delay coordinate 

t 0 turbulent timescale 

l turbulent lengthscale 

k turbulent kinetic energy 

s turbulent dissipation rate 

/ip two-point, fourth-order space-time correlation of velocity 

Rij two-point, second-order, spatial correlation of velocity 

h temporal correlation function 

H Fourier transform of temporal correlation function h 

N non-compactness factor 

c t calibration constant associated with lengthscale 

c z calibration constant associated with timescale 

A m calibration constant associated with conversion of turbulent kinetic energy to acoustic energy 
a calibration constant to remove singularity of Doppler factor 
F master spectrum, spectral function at d = 90° 

G , G Green’s functions 
()c () convecting 

()oo () at far-field 

()’ fluctuating portion of () 

() s () evaluated at source location 

() t () evaluated at transverse location 

[] m [] in moving frame 

( ) time average of () 
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2. Introduction 


There is an ever-increasing demand on the aeroacoustics community to achieve new levels of 
accuracy in noise prediction in order to meet the stringent regulations set for aircraft noise reduction. Jet 
noise is of particular interest because of its dominant contribution to the noise radiated from aircraft 
engines operating at high thrust conditions during takeoff. As part of an ongoing research program at the 
NASA Glenn Research Center, the JeNo code is being developed to provide the state-of-the-art in jet 
noise source modeling and propagation. Since jet mixing-noise is broadband, with frequencies that cover 
a full three-octave range, a successful source model needs to bring in enough physics to predict a 
reasonably acceptable spectral density at various operating conditions. 

This document, Part 1: JeNo Technical Manual, details the theoretical development of JeNo version 
1.0. It discusses the jet physics being captured by the code, the mathematical derivation of the equations, 
the assumptions and models which are captured in the computer code. This document is for the user who 
wishes to understand the limitations of JeNo with either an eye to understanding the output or to 
improving its shortcomings. A separate document, Part 2: JeNo Users Manual, covers the operation of 
the code from the computer user’s point of view, including code architecture, input and output formats, 
and examples. 


2.1 Jet Noise Source Physics from Observations 

Experimental observations of jet noise spectra and directivity for a wide range of shock- free operating 
conditions show that a bi-modal description best fits a typical spectrum. This has prompted many to 
propose a two-source generation mechanism. In the broadside angles, jet noise is usually attributed to 
small-scale turbulence. Prediction schemes such as MGBK(ref. 1) or Tam and Auriault (ref. 2) scale the 
sound spectral density to 7/2 power of turbulence kinetic energy and use a Reynolds-averaged Navier- 
Stokes (RANS) solution to estimate the source strength and its spectral shape. The general shape of the 
spectrum and its roll-off at the high- and low ends of its frequency range depends on the proposed source 
model as well as the accuracy of the propagation filter, i.e., the Green’s Function (GF). The non-causal 
GF solution of interest may be obtained from the inhomogeneous Lilley’s equation or directly from 
linearized Euler equations using numerical methods or some form of analytical approximation. This is the 
source being addressed in version 1.0 of the JeNo code. 

At aft angles, the non-causal GF solution predicts a strong increase in sound, but not as strong as 
observed in jet experiments at supersonic Mach numbers. A second mechanism for sound generation is 
commonly proposed, being attributed to the large-scale turbulent structures, which are considered as a 
superposition of instability modes of the mean flow. It has been shown that when the base flow is 
properly represented as a non-parallel jet, these waves grow in amplitude and then decay farther 
downstream. They become increasingly more efficient with jet speed, and at small angles from the 
downstream axis. While this mechanism is often called a second source, a recent study (ref. 3) describes 
instability waves as a second ‘conduit’ that helps carry the sound from the source region to the far field 
rather than acting as a separate source. In this view, a two-component spectrum may be thought of as the 
contributions from two completely different Green’s functions acting on the same source. This second 
‘conduit’, expressed by a causal GF, is not addressed in version 1.0 of the JeNo code. 

To summarize, version 1.0 of JeNo is applicable to cold, round jets of subsonic velocity. It correctly 
predicts sound for all angles at such conditions. It does not correctly predict hot jets, jets of strong 
asymmetry, or jets with supersonic (relative to ambient speed of sound) exit velocities. 
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2.2 Outline of JeNo Theoretical Development 


The non-causal GF described in the present paper is derived from Lilley’s equation. A good 
approximation to the base flow that is useful for jets and adopted here is a unidirectional transversely 
sheared mean flow or a locally parallel flow. From a practical point of view, this flow provides added 
numerical efficiency with minimal degradation in accuracy. Deviations from a true spreading jet are 
usually limited to a small angle range close to the jet axis and near the zone of silence where the 
contribution from a causal GF is likely to dominate. 

Upon linearizing the Navier-Stokes equations about the base flow (for an ideal gas), the nonlinear 
terms are moved to the right-hand side to obtain the inhomogeneous Pridmore-Brown equation 
(ref. 4). The right-hand side terms, all second-order in fluctuating variables, are now identified as the 
equivalent sources. Fluid viscosity and heat conduction are considered to be relatively unimportant in 
noise generation and are neglected compared to other source terms. Two of the three remaining sources 
are quadrupole in nature and are referred to as the self- and shear-noise source terms. The third source has 
a dipole character and is associated with temperature fluctuations. 

Solving the equations for the far-field observer then consists of two tasks: modeling the source terms 
and finding the solution to the Pridmore-Brown equation. The first is done by assuming functional forms 
for the different terms and applying proper scaling factors from turbulence models, usually obtained from 
CFD. The second is done by constructing a Green’s function (GF) which satisfies the propagation 
equation for a point source embedded in the flow, which is also usually known from CFD. This GF then 
acts as a filter to describe the relationship between the sources and the observer. 

A model-based approach usually relies on averaged equations of motion. Statistical properties of 
noise sources, as described by a two-point space-time correlation, are entirely modeled. A physics-based 
modeling of sources can be improved as our understanding of the underlying generation mechanism and 
turbulence details advances. The current source models represent current understanding of the turbulence 
correlation tensor and will undoubtedly be improved upon in future versions. It is important that the 
source model be kept independent of propagation effects so that improvements in the source model not be 
masked by improper fixes to the propagation portion of the code. 

In addition to the source description, the propagation filter needs to be accurate enough to capture 
refraction and shielding of sound for a wide range of frequencies, angles, and operating conditions. A 
review of some commonly used high-frequency asymptotic forms of the GF following the usual WKB 
methods and comparison with the numerical solution (ref. 5) demonstrates that these solutions generally 
remain accurate down to a Strouhal number of 0.5. However, the extension of the WKB solution to the 
more general jets with multiple turning points presents additional challenges and may require a numerical 
calculation of the near-field solution in order to satisfy the matching conditions. Indeed, the initial version 
of JeNo differs from its ancestor, the MGBK code, primarily in its use of the more rigorous GF solution, 
obtained using adjoint methods. 


2.3 Outline of Document 

The objective of this paper is to present the details of the source as well as the GF in a model-based 
prediction approach. The analytical solution to the intensity calculations using an exact GF is given in 
section 3. Section 4 shows the modeling of the two-point velocity correlations using an exponential form 
to describe both spatial and temporal functions. The general features of the 90° spectrum are addressed in 
section 5. A parametric study of the GF using a RANS-based mean flow is discussed in section 6. In 
particular, it is demonstrated that the mean flow could have an amplifying effect on radiated sound at 
certain observer angles and at certain source locations. More importantly, it is argued that regardless of 
the source definition used in the convolution integral, a non-causal GF, due to its sizable zone of silence, 
is incapable of capturing the observed directivity of high-speed jets at their peak radiation angle. Sample 
noise predictions are presented in section 7 for a number of unheated jets in the Mach number range 0.5 1 
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through 1.80, and compared with data recently collected at the Small Hot Jet Acoustic Rig (SHJAR) at 
the NASA Glenn Research Center. Section 8 studies a Mach 0.90 heated jet at a stagnation temperature 
ratio of 3.14. The hot jet exercise considers two important scaling rules related to the GF, and proposes a 
number of turbulence-related issues that need to be examined carefully prior to the addition of a heat- 
related source. A summary of assumptions adopted in the present analysis will be reviewed in section 9. 
Requirements for developing a comprehensive physics-based jet noise prediction algorithm will also be 
discussed. 


3. The Governing Equations 

Lilley’s third-order wave equation (ref. 6) may be linearized about a unidirectional transversely 
sheared mean flow 

Po = constant, U / = S a U (x 2 x 3 ), T 0 = T 0 (x 2 x 3 ) . 


and rearranged to a form commonly referred to as the Pridmore-Brown (ref. 4) equation 


LU = — 
Dt 


V 


p-n 

Df 


dx 


■(c* an 


dx 


+ 2 c 


j ) 


du d 2 n 

dxj dxfixj 


= r 


(i) 


Here T is the source, x = (x l ,x 2 ,x 2 ) are the Cartesian coordinates and x\ is in the stream-wise 
direction, p is pressure and p a is the mean pressure, (i.e., p = p o + p' ), FI = (1 If) In (p/p 0 ), 
c 2 = ;/.R T„ is the mean sound speed, t is time, and the convective derivative operator is 

D _ d | d 
Dt dt dx j 


When pressure fluctuations p' are small relative to mean pressure p a , the dependent variable is 
approximated as II = p'/(/ p o ) . Goldstein (refs. 7 and 8) carried out a second-order expansion of 
equation (1) and introduced a new dependent variable n'= II + II 2 / 2 to show that 


Dt dx i dXj dx l 


(2) 


and 




TtM +(c2) ,sn 

dxj dx t 


(3) 


where «, = v,—SnU denotes the velocity fluctuations. The source term (3) may be approximated as 
Ci = d(uiUj)/dxj once fluctuations in sound speed (c 2 )'= /)R(T-T o ) are neglected (a valid approximation 
when jet is isothermal). With the dependent variable approximated as n'=Yl 
(see ref. 9) for the equation governing FI 2 ), equation (2) becomes 
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p 


( 4 ) 


TT1 ~ D 5 2 ( m , m 7 ) ^ dU d 2 (u ,Uj) 

XjI 1 = Z ^ XX 

Dt dXjdXj dx i dx\dXj y p 0 

An alternate definition for the dependent variable is 7r'=(p / p o )' : — 1. As shown in (ref. 7), with 
this new parameter, equation (2) becomes an exact rearrangement of the Euler equations once the source 
is defined as 


Ci=^-0-+^)u t Uj+(r-i) h'^- 

OXj ox t 

This latter source reduces to that defined in equation (3) when enthalpy fluctuation h' is written as 
h '=(/ — l)(c 2 )' and if n' is considered small relative to unity. 

The equivalent sources of sound in a Lilley-type analogy are second-order in fluctuating variables. It 
has been pointed out repeatedly in the literature that the main difference between Lighthill and Lilley’s 
formulation of the acoustic analogy is a separation of the propagation effect from the source in the latter 
formulation, which presumably places less demand on the details of the source modeling. The GF 
calculations described in section 6 clearly show the significance of the Lilley operator in isolating the 
propagation effect. 

There have been numerous efforts in the past half-century to experimentally measure the actual 
sources of sound in a turbulent flow. Measurement techniques such as phased-array can measure the 
relative source strength and its frequency content at various regions in a jet. Particle Image Velocimetry 
(PIV) methods have been instrumental in turbulence and correlation measurements. When combined, the 
two techniques may establish a link between the turbulent flow and the far-field sound it produced. This 
rather laborious and cumbersome process is used in concept studies and parameter selection. 

The analysis presented in this paper is based on equation (4). The 2 nd order terms on the right side of 
the above equation are referred to as self- and shear-noise source terms, respectively. 

3.1 Green’s Function 

A stationary point source with frequency co and location x s (superscript s denotes a source location) is 
considered in defining the GF 


L(&e‘ a, ) = cle w, S(x-x) 


(5) 


Using an adjoint method (ref. 10), the far-field GF is given as 


G(x,x ! ,co) 


1 gik(R-x{ cos G) 

\ncoR 


m = 0 


( 6 ) 


where ( R,0,<p ) are the observer spherical coordinates, with radius R measured from the jet exit 
centerline, polar angle 9 measured from the stream- wise jet axis, and azimuthal angle cp is measured in a 
span-wise plane. Wave number is defined as k = col c* and r = R sin 9 Function f„(r, k, 9) is a solution to 
the second-order compressible Rayleigh operator subject to appropriate matching conditions at the jet 
boundary (app. A). The linear operator L also supports instability waves, which for a jet are the well- 
known Kelvin-Helmholtz instabilities. Reference 1 1 argues that these waves are suppressed if the 
governing equations are solved in a frequency domain and if a time-harmonic response is assumed. 
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For a source typeD/ Dt^c^e ' M 8(x — X s )], representative of the self-noise term in equation (4), the 

GF is denoted asG, and is related to 6 as 6 = —{io) + Ud I dx[)%, , which upon using 
equation (6) becomes 


G = -ioj(\ - M s cos#)G 


(7) 


Flere Af = U{f)/c m is the acoustic Mach number at the source location. 

The GF of interest in jet noise is associated with a moving singularity with source frequency af and 
convection velocity i U c . For self-noise we write 


D 


L(Ge^ = ^{cle-^S( Xl -U c t)S(x t -x;)} 


( 8 ) 


Subscript t denotes a transverse location. The above GF is derived from equation (7) and a 
convolution integral as 

G(x,x s ,co)e uot = — [ dco \ [ G(x, y, co)e^ icot e l(co ~ co)T 8{y t - x^)8(y t -U c r)drdv (9) 

2 7T JJJ 


co y t 


Integrating with respect to r and y, yields 

G(x,x s ,co)e Uo1 = dcoj&(x;x;,y l ;(o)e- M e- i(6>s - ,aM dy 1 


(10) 


c co y x 


Using G from equations (7) and (6) into equation (10) and invoking the parallel flow assumption; the 
integration with respect to y t yields a delta function 

f = 2nS[ _ _^_ (1 _ ^ CQS ^ + ^ (1 

i U c U c 


Where M c = U c /c x is the convective Mach number. 

Placing equation (11) in equation (10), the integration with respect to co is readily carried out 


G{x,x s 


• /i n /fS /i\ ico t ico / Cqq si 

^ (1 M COS 0) ^\-M c cos0 g\-M c o,osO \ 1 r ss ® ' ^oo 

AxR M c cos6) m l-M f cos#’ 


6) cos m{cp - cp s ) 

( 12 ) 


Upon relating source and observer frequencies through the Doppler factor 


co = 


co 


1- cos# 


(13) 
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the GF to equation (8) becomes 


G(x,x s ,co) 


-i (l-M'cos#) 
4 nR (1- M c cos6) 


e ' kR Z fm( rS ’ k ’ &) C0S m (<P ~ <P* ) 

m 


(14) 


When the source is moving, the Doppler factor becomes singular in the direction 6= cos ’(1 /M c ). By 
considering the source as non-compact (app. B), a modified Doppler factor is defined that should be used 
near this singularity. 

Equation (14) indicates that the adjoint equation should be solved at the observer frequency even 
though the source is moving. In addition, it illustrates that for a convecting type source, the dependence of 
the GF on the axial source location is purely implicit, i.e., through the flow definition at the jet slice. Once 
the adjoint equation is solved for given values of observer angle 9, Strouhal number St, and mode number 
m, function f m should be known at any arbitrary source locations r on the jet slice. Strouhal number is 
defined using exit diameter and jet exit velocity St = coDKlnU j). 

3.2 Multi-Pole Sources 

The far-field spectral density may be expressed as the integration of the sound spectral density per 
unit volume of turbulence with respect to the source volume y (app. C). 


p 2 (x ,y,a>) = J J G\x,y-} / 2,co)G(x,y + } /2,a>)R(y,},r)e U0T dTd} (15) 

| -oo 

Herei?( v,£ r) denotes a two-point fourth-order space-time correlation between source points 

y t = y - q S 2 and y 2 = y + g/2 separated by time r and space £ . In a compact eddy approximation, the 
variation of the GF within the source region is usually neglected and the product GG* is approximated as 
the square of the magnitude of the GF at the center of the correlation y, and taken out of the above 
integral. The phase variation in the GF, however, could become a factor at high frequency when the 
correlation length scale exceeds the acoustic wavelength. The complex function f m appearing in the series 
summation of equation (14) contains the phase information of the GF. To simplify the analysis, a form of 
the phase factor applicable in the high-frequency limit is used here (app. D) 

+00 

p 2 (x,y,a>) = \G(x,y,co)\ 2 j fe ,k4 R{y,^,r)e ,0>T dT(f^ (16) 

l-CO 

The above double integral, referred to as the wave-number frequency spectrum, may be expressed 
either in a fixed or a moving frame. In a convecting frame, iq m = %—iU c T, and 
R(y, £ T) = R m (y,i m , 0, hence 

+00 +00 

j e'-'drj R(V,1 r K‘W|= j e“’’dr\ (17) 

— 0 ° g —00 
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A spectral coefficient for a two-point fourth-order correlation of stress component u ; u, at points y 
and y is defined as (prime designates separation in space and time) 

+co 

I ijk ,(y, co) = | J (u^jXu'bU,) e ,k £ e l0,r circle ( 18 a) 

| -00 

Upon recognizing that —k . % + (OX = —k . g m + of X , the above integral is written in a moving frame 

+00 

hjkAy^) = J J [(^UjXu.u,) l n e~ ,k4m e U0 ' T dxd} m ( 18 b) 

Subscript m under the square bracket in equation (18b) denotes a moving frame two-point correlation, 
which is related to the nozzle-fixed correlation by an exchange of £\—U c x in the latter with (£l) m . The over- 
bar points to a time-averaged quantity. Presence of double derivatives operating on the source terms of 
equation (4) indicates that the appropriate GF should ultimately be calculated for a singularity type 

Dt ox fix . 

1 J 

The analysis may be carried out as described in the convolution integral leading to equation (9). 
Flowever, the source spatial derivatives are now moved to the GF prior to y integration. For example, an 

application of 8" / dy l dy l to G(x, v; co) in the integrand of equation (9) generates (—ik cos#)” as seen 
from equation (6). Subsequently, (kcosdf multiplies sound intensity due to source 
correlation (m 1 m 1 )(m 1 w 1 ) . A derivative of G with respect to the span-wise source coordinates y 2 or y 3 , acts 

on f„(r\k, O)cosm( (p-(p). The algebra is quite lengthy but relatively straightforward. It may be shown that 
to the first order of approximation, quadrupoles of self noise, once integrated azimuthally for a ring 
volume within an axisymmetric jet, relate to the axial component 7 ini . For an isotropic turbulence the 
result simplifies to 


(\-M s cos 6)\p s ) 2 k A I in 


( 20 ) 


This, in view of equation (14), leads to the following expression for the sound spectral density per 
unit ring volume at radius r' 


p se if(x,y,OJ) = 


1 (l-ATcos#) 6 ^ 


( 4 txR) 2 ( 1 - M c cos0) 


P_ 

VPoo) 




(21) 


The above result is integrated over jet volume elements that radiate directly to the far field. Noise 
from sources internal to the jet requires special treatment due to reflection from solid boundaries. An 
approximate high-frequency solution for sound emission due to internal mixing is given in 
reference 12. 

It will be shown in section 7 that equation (14) scales like /? 1/2 with jet density, hence sound spectral 
density in equation (21) scales with p\ 
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4. Source Model 


The axial correlation coefficient needed in equation (21) is obtained from equation (18) when all 
indices are set equal to one. Using the usual quasi-normal approximation for the joint probability 
distribution of turbulence (refs. 13 and 14), a fourth-order correlation is written as a product of second- 
order tensors. 


u i u j u k u i = (u i u J )(u' k u ( ) + (u i u' k )(u j u { ) + (u i u l )(u j u k ) 
A nozzle-fixed second-order correlation is written as 


UiUj (i,r) = Rj(}- i U c T)h( t) 


(22a) 


(22b) 


Appropriate modeling of a two-point velocity correlation is a crucial step in a physics-based 
prediction approach. As shown in reference 1 the best representation of the experimental measurements 
(refs. 15, 16, and 17) are obtained when exponential functions are selected to represent both spatial and 
temporal dependencies of the correlation (also see app. E) 

AU = e-^ le , h(T) = e+ ,r - 1 (23) 


In equation (23), A I I is measured in a moving frame, and is equivalent to 

~ ^ J) + + + A, ] if used in a fixed frame, and f is the scalar function that appears in the isotropic 
from of Rij. 

In an isotropic turbulence, various components of a two-point correlation are obtained from placing 
equation (23) in Batchelor’s isotropic turbulence model ref. 14)-for example, the fixed-frame axial 
component of the correlation is written in the following non-separable form 


n 






4=[(&-u e rf + Z + g] 1 


Parameters ^ and to denote the length- and time-scales of the correlation 

i = C ( K V1 la, r o = c t k / £ (24) 

Turbulence kinetic energy k and its dissipation rate s are usually provided through a RANS type 
(c cl 

solution and v f ’ 1 represent a pair of calibration constants. 

Using a quasi-normal approximation, and substituting equation (22b) in equation (18b) shows that 

4m 


+00 

H(co s )= jh 2 (T)e“° T dz 


-00 


(26) 
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Upon representing R n (£ m ) according to the homogeneous isotropic model of Batchelor (ref. 14) 
with / as defined in equation (23), it is shown (ref. 1) that 


Am (y**® ) 


— Y (uy£ 3 H(o) s )N(k £ ) 

5 71 


(27) 


1 1 ere N(kt) denotes a non-compactness factor (app. F), which appears in the analysis due to the 
phase variation of the GF. As was discussed in reference 1, N(k£) remains equal to 1.0 when 
0 <ki< 2n, and decays rapidly as the eddy length-scale exceeds the wave-length of acoustic 
disturbances (i.e., ki > In'). 

Carrying out the integration in equation (26) in conjunction with equation (23) results in 


H(co s ) 


\ + {co s r 0 / 2) 2 


(28) 


Upon using equations (24) and (28) in equation (27) and replacing u x with 2 a'/3 gives 


16 ,c 


hC’/ 2 


Am (y,a s ) = jrr(—) « 
45;r c 


\+{G) s T o ny 


- N(k £ ) 


(29) 


Equation (29) shows that at constant r„ sound intensity scales with a 7 ' 2 . A third calibration constant 
A m should multiply equation (29) to account for the fraction of turbulence kinetic energy converted to 
sound. This factor is combined with other constants preceding a 7 2 


Am (y, 0 ) s ) = A m K 


7/2 


\ + (co s r 0 / 2)' 


-N(k£) 


(30) 


Equation (21) in conjunction with equation (30) is integrated over the jet volume (all rings within 
each slice) to calculate the far- field sound due to self-noise sources. 

The appropriate GF for shear-noise source term is similar to equation (14), but divided by 
oi(\ -M s cos 0). Once the variations in dU/dr within the source correlation volume are neglected, the shear- 
noise spectral density due to a unit ring volume may be approximated as 


( 


2 dU / dr ^2 




co 


(^-cos 2 0)1 


(4 nRf (1- M c cos0) 2 '2 

<(^Am)^Z|.A,(^^^)| 2 - 


cos 2 9 + 


(1- M s cos6*) 

c-l cl 




(31) 


This expression remains small at and near mid angles; however it generates significant directivity 
away from 90° as opposed to the self-noise spectrum, which is more omni-directional. 

The coupling between self- and shear-noise terms is neglected in the following predictions. 
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5. Spectrum Function 


The spectral shape at 90°, also referred to as the master spectrum, is derived from the product 
k 4 H(a>)N(kl), as seen from equations (21) and (30). It takes the following form once turbulence 
correlation is modeled according to equation (23) 

F(a>z„) - (tt,r ‘. )4 2 N(kf) (32) 

1 + (gw 0 /2) 

In the absence of factor N(kl), function F(oot 0 ) becomes infinite like of as < 0 —>oo . On the other hand, 
with increasing co the phase variation in the GF becomes increasingly significant as the eddy correlation 
length exceeds the acoustic wavelength. The spectral function F(cot 0 ) decays at large frequencies as of 3 
once factor N(k£) is included (fig. 1). This practically amounts to limited self-cancellation of high wave- 
number components of noise within the correlation. Also shown in this figure is the 90° spectrum 
function as calculated with temporal function h(r) = exp[-^(<j/ 2) 2 +(t / r o ) 2 ]. The small constant cr 
would assure high-frequency decay even with factor N(kl) selected as 1.0; however the spectrum 
becomes narrower. A Gaussian temporal function cxp(-r7zy ) generates a relatively narrow spectral 
range due to its steep roll-off at high frequency and may not be a suitable model for jet noise prediction. 


6. Numerical Evaluation of the Green’s Function 

In this section numerical calculation of the GF is presented for a range of parameters of common 
interest in jet noise prediction. A careful study of this non-causal GF and its association with the source 
location, frequency, and jet speed provides valuable insight into its contribution to a total solution. 

In a locally parallel flow, jet profiles (i.e., mean axial velocity and temperature) as defined at each 
stream-wise location X\, are utilized in solving the propagation equation (Al), subject to the boundary 
conditions defined in appendix A. The GF is calculated from the series summation equation (14)-and 
subsequently normalized and integrated azimuthally to define ring source directivity 

D 2 (x,x s ,co)= J|G(jc,jc J ,fiO| d( P S (33) 

-71 

Both stationary and convecting type sources are considered. Source convection velocity is defined as 
a weighted average of local and exit velocities as supported by measurements (ref. 18) 

U c (r,x l ) = 0.50(7 (r,X\) + 0.25(7,. (34) 


Alternative expressions (ref. 19) for U c based on local, exit, and ambient jet velocities are available 
and could readily be tested. When the denominator of equation (14) becomes singular, it is replaced with 
a modified Doppler factor (app. B). 

Mean flow predictions (RANS) were generated for a Mach 0.90 unheated jet using the Wind flow 
solver (ref. 20) available at the NASA Glenn Research Center. Source location y is identified through the 

ring radius r and its distance from the jet cxitx,' . 

Directivity D for a stationary ring source at one diameter from jet the exit, and for selected values of 
St from 0.10 to 5.0 is shown in figures 2(a) through (d). It is observed that sources in the proximity of the 
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jet centerline produce a larger directivity. The zone of silence grows noticeably in size with increasing 
frequency as expected. As seen in figures 2(c) and (d), the zone of silence may grow as large as 60° 
(relative to the jet axis) for the more energetic segments of the jet, i.e., the jet mixing layer. Figures 3(a) 
through (d) use the flow definition at x ,' =7.5/1 to highlight the significance of the jet profile on sound 
refraction. It is noted that source strength should ultimately multiply the GF when calculating the actual 
sound spectrum. 

Source convection is known to amplify the radiated sound in the mean flow direction. This is 
illustrated in figures 4 and 5, which should be compared with the figures 2(b) and 3(b) for similar 
stationary sources. 


7. Acoustic Results 

A total of seven unheated cases (stagnation temperature ratio of 1.0) are shown in table 1, with the 
aerodynamic Mach number extending from 0.51 to 1.80. All nozzles have a 2 in. exit diameter. The first 
two set points are subsonic-and utilize an Acoustic Reference Nozzle (ARN2). The supersonic set points 
use Convergent-Divergent (CD) nozzles designed to provide perfectly expanded exit conditions using a 
Method of Characteristic approach. The fully expanded (design) Mach numbers are 1.185, 1.40, 1.50, 
1.66, and 1.80. Jet mixing noise was predicted on an arc RJDj = 100, and compared with SFIJAR data 
recently acquired at the NASA Glenn Research Center. Atmospheric attenuation has been removed from 
all measurements in order to make a lossless comparison with predictions. Although special care has been 
taken to test each CD nozzle at its design pressure ratio, shock-associated noise and screech-related spikes 
are noticeable in the measurements. RANS solutions were generated using the in-house Wind code with a 
standard k-e turbulence model. Upstream nozzle conditions were specified in terms of plenum 
temperature ratio T r and the pressure ratio. 


TABLE 1.— UNHEATED JETS (T R = 1.0) 


Case 

Mach No. 

Uj/C K 

SP03 

0.51 

0.50 

SP07 

0.98 

0.90 

Ml. 18 

1.185 

1.047 

Ml. 40 

1.40 

1.186 

Ml. 50 

1.50 

1.245 

Ml. 66 

1.66 

1.33 

Ml. 80 

1.80 

1.40 


TABLE 2.— HEATED JET (1^ = 3.14) 


Case 

Mach No. 

Uj/c x 

SP49 

0.90 

1.48 


An additional heated (nearly sonic) case, SP49, will also be considered to highlight the shortcomings 
in our hot jet predictions and to suggest possible solutions. 

Calibration constants c T ,c f (see eq. 24) and A m (see eq. 30) are determined for the best fit with respect 
Mach 0.51 jet. The fidelity of the solution is then examined for the range of Mach numbers. The argument 
of the non-compactness factor ki = {2nSt\t / D ; )(Uj / c t ) would become self-similar across a wide 
range of Mach number when divided by the local acoustic Mach number U(r*)lc x , as both factors t / D j 
and Uj/U(r ) are now self-similar. This modified argument is deployed in noise predictions. 

Figure 6 shows the spectra for a Mach 0.51 cold jet at selected inlet angles. The self-noise component 
exhibits an omni-directional appearance as a result of the factor (1-AT cost?) 6 in the numerator of 
equation (21). The shift in the spectral peak with angle is primarily due to the mean flow effects and is 
captured by the GF. Figure 7(a) shows similar comparisons in a Mach 0.98 jet. The agreement is 
generally good, both in level as well as the spectral peak location. Near the dominant radiation angle of 
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the Mach 0.98 jet, however, differences start to emerge. The data point to a distinct spectral peak near the 
150° inlet angle and at the Strouhal number of 0.20. It is suggested here that this peak may be influenced 
by the Kelvin-Helmhotz instabilities, which manifest themselves at high subsonic and supersonic Mach 
numbers. These waves dominate the distinct directivity of high-speed jet noise at shallow angles, and 
their peak radiation frequency is sometimes scaled with Helmholtz number (He = fDJcf). 

Similar comparisons are also shown in figure 7(b) for the five supersonic jets. The presence of shock 
noise in the spectra is quite visible. These predictions exhibit increasing deterioration at shallow angles. 
The 150° spectra peak at a lower frequency relative to measurements and show a deficit in noise level. It 
is conceivable that, in addition to the instability -related noise, jet spread could also play a role at high 
speeds and near the zone of relative silence. A parallel flow model tends to exaggerate the High 
Frequency (HF) refraction (refs. 21 to 23). The actual decay rate into the zone of silence is relatively less 
steep when the HF-GF is calculated in a spreading jet. The effect of jet spread on Low Frequency (LF) 
noise is less clear and yet to be examined. For long waves, the locally parallel flow assumption is more 
likely to deteriorate when the mean flow is not slowly varying on a wavelength scale. On the other hand, 
it is also argued that refraction is relatively weak at low frequency. 


8. Hot Jets 

When the acoustic Mach number exceeds the sonic point, the factor (1 —Af cost?) appearing in the 
denominator in the compressible Rayleigh equation becomes singular. The numerical integration should 
now continue along an appropriate path in the complex plane, and around the singularity 
(ref. 10). A better illustration of the GF that shows its scaling with respect to square root of density is 
defined in equation (35). The factor (1 -Af cost?) is also multiplied by the GF to demonstrate the transition 
into the zone of silence at cost? = 1 /Af. 

A(x,x s ,co) = (\- M s cos 0)(p s / p^y 05 D(x,x s ,co) (35) 

Factor D 2 was defined in equation (33). As noted above, the appearance of the density factor in 
equation (35) solely serves to show that the GF scales with the square root of density at moderate to high 
frequency. This is examined by observing the behavior of the GF at 90°. 

Figures 8(a) and (b) show the GF for a stationary ring source in a Mach 0.90 hot jet, and at a 
temperature ratio of T r = 2.78. At this temperature the acoustic Mach number is supersonic 
(i.e., U/c x = 1.45). Inspection of the GF at 9= 90° and at St= 0.10 shows some dispersion with the source 
location (fig. 8(a)). Figure 8(b) shows that the above density scaling is successful at St = 1.0. More 
importantly, this figure demonstrates that in the mid- to high frequency range, the GF equation (14) scales 
like (l-Af cost?) -1 outside its zone of silence 

It is also evident that compared to unheated subsonic jets discussed earlier, the zone of silence has 
grown in size irrespective of the frequency (see figs. 3(a) and 8(a)). The implications are that the above 
GF is not going to provide the required noise directivity near the peak radiation angle. Figures 9(a) and 
(b) show similar results for a convecting ring source. Figure 10 compares the 90° spectrum in a 
Mach 0.801, T r =3.0 hot jet as calculated with an exact GF as well as a density-scaled GF (i.e., GF is 
simply set equal to p m /p 12 ). The difference is fairly small and limited to the low frequency range. 

The spectrum at R/Dj = 100 at SP49 (table 2) is also compared with the SHJAR data as seen in 
figure 1 1 . Consider the 90° spectrum, which exhibits a deficit of 3 to 8 dB relative to measurements. This 
spectrum is entirely defined through equation (21), which according to the preceding arguments scales the 
sound spectral density as/? 3 , in addition to the usual k 7 ' 2 scaling. Therefore, the mean density reductions 
caused by heat addition need to be compensated for by other resources in order to reach the required noise 
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level. The choices available are: (a) additional noise sources due to heat; (b) changes in turbulence and/or 
its scales; or some combination of (a) and (b). 

Lighthill’s formulation of the acoustic analogy shows that in the absence of the viscous stress tensor, 
the source consists of an additional term (p'-c^p'). When entropy fluctuations are neglected, this 
becomes ( c 2 -c 2 x )p ' which is not necessarily small compared to terms usually kept in the equation. 
Unfortunately/?' now appears on both sides of the equation. Other forms of the analogy have tried, for 
example, new definitions for the dependent variable, in an effort to restrain the impact of those 
inconvenient terms that are usually neglected. It is argued (ref. 7) that the second term on the right hand 
side of equation (3) is of the order of ( U / c) 2 (u 2 / /) relative to the first term that scales like ( uli ). In the 
Mach 0.90 example seen in figure 11, the second term has a magnitude of 0.81 compared to the first, 
which makes it relatively significant. Modeling of this term, however, poses new difficulties due to the 
appearance of p' in the source. 

In the latest formulation of the acoustic analogy, reference 24 selects a new dependent variable, 
referred to as the generalized pressure, and writes the third-order wave equation in a form that does not 
include the dependent variable as part of the source. 

Tam et al. (ref. 25) modified their ks turbulence model in order to improve their fine-scale predictions 
for hot jets. They proposed a new addition to the turbulent viscosity, as related to the mean density 
gradient, which accounts for enhanced mixing and increased growth of the shear layer when jets are hot. 
More recently (ref. 26), they modified their two-point space-time correlation model to increase its decay 
rate for hot jets. This new model, which takes advantage of four additional empirical parameters, provides 
improved prediction for hot jets at mid-angles. 

Shown in figures 12(a) through (c) are the predicted turbulent kinetic energy for the 2 in. convergent 
nozzle at set points SP03, SP07, and SP49 (see tables 1 and 2), using a modification to the baseline k-s 
model in the Wind code that incorporates variable diffusion coefficients (ref. 27) in the transport 
equations for k and s. The objective of this modification was to enable more accurate calculation of the 
jet potential core region, and in particular, the faster transport of turbulent kinetic energy to the jet 
centerline observed in experiments, which has not been captured well by standard k-s models. Using this 
modified k-s model, favorable agreement with the PIV measurements (ref. 16) are observed in the 
turbulence levels, as well as the shortening of the core length due to the heat-which is the product of a 
faster jet spread and improved mixing. Consequently, it is quite compelling to attribute the noise deficit 
observed in the predictions of figure 11 to heat-related sources, and/or changes in the turbulence scales, 
which in turn, may influence the decay rate of a two-point correlation. 

Near the downstream jet axis, the instability-related noise quickly overwhelms the hot jet solution. 
The problem is also compounded due to a gradual deterioration in the parallel flow model at higher jet 
speeds and at small angles, as discussed in the latter part of section 7. 


9. Assumptions and Future Work 

This paper details some important features of the source and the propagation Green’s function in a 
model-based jet noise prediction methodology. A number of hypotheses are made in calculating the far- 
field sound: 

A careful examination of the propagation filter shows that in a locally parallel flow the zone of 
silence grows increasingly large with jet velocity. It could encompass a sizeable region near the 
downstream jet axis as the acoustic Mach number becomes supersonic. A parallel flow model is known to 
increasingly exaggerate the high frequency refraction at higher jet speeds (refs. 21 and 23). When the GF 
is calculated in a truly spreading jet, the actual decay rate into the cone of relative silence is less steep. 
Low frequency noise may also be influenced by the jet spread-although the effect is less clear and needs 
to be determined. In general, a parallel flow model is expected to provide reliable sideline predictions, 
and to become effective in a broader angle range at lower jet speeds. 
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TABLE 3.— LIST OF JENO ASSUMPTIONS 



Assumption 

Reference 

1 

Turbulence is assumed as quasi-nonnal in order to reduce the order of the fourth order 
correlations 

eq. 22(a) 

2 

Fluid viscosity and heat conduction are not considered to be important in sound 
generation. 

eq. 1 

3 

Turbulence is considered as isotropic 

eq. FI; also ref. 14 

4 

Exponential spatial and temporal functions are utilized to model a two-point correlation 
function 

eq. 23 

5 

A non-compact source is considered 

eq. 27 

6 

The non-causal GF is calculated for a locally parallel flow 

eq. 1 

7 

The second component of the GF, believed to be associated with the instability waves, is 
not addressed 

see ref. 3 

8 

The heat-related sources were not modeled 

see equation after eq. 4 

9 

Two constants are utilized in connection with the time-and length scales. 

eq. 24 


Third constant relates turbulent kinetic energy conversion to sound. 

eq. 30 


Fourth constant removes singular behavior of the Doppler factor. 

eq. B6 


The second component of the Green’s function (i.e., the causal GF) is instability related, and 
dominates the spectrum at shallow angles to the jet axis. This component needs to be addressed as a high 
priority in supersonic jets. Some preliminary results calculated for a simplified source in a model jet 
(ref. 3) show encouraging qualitative agreement with data. 

It was also shown that a non-compact source model predicts a proper high-frequency roll-off. As the 
wavelength of the sound generated within a turbulent eddy becomes shorter than the correlation length- 
scale, a phase cancellation is initiated that shapes the spectral roll-off at high frequency. It is not 
surprising that this fall-off should be sensitive to the spatial form of the correlation in relationship to the 
wavelength of the sound it generates. Various spatial functions such as Gaussian versus exponential 
become distinguishable from one another when the source is considered as non-compact. 

Although the formulation and validation cases presented here are limited to unheated jets, we took the 
liberty to evaluate a sample hot case to demonstrate the shortfalls and possible solutions. The deficit in the 
predicted noise level was attributed to any number of the following: 

• Calibration of the usual two-equation turbulence model. 

• Possible additional sources such as enthalpy fluctuations. 

• Decay of a two-point space-time correlation as affected by heat. 

• The effect of heat on turbulence length- and time-scales. 

Experimental studies are currently underway at the NASA Glenn Research Center to address items 
(c) and (d) above. In addition, turbulence model development work is ongoing in the Wind code with the 
addition of linear two-equation and explicit algebraic stress models specifically calibrated for jet flows. 
Flowever, the favorable comparisons of P1V measurements with RANS calculations using a modified k-s 
model presented here indicate that the RANS input is not the sole source of disagreement between 
calculated and experimentally measured jet noise. It is hoped that these efforts may shed light on this 
somewhat elusive relationship between the turbulence in hot jets and its far- field noise. 

To reiterate, version 1.0 of JeNo is applicable to cold, round jets of subsonic velocity. It correctly 
predicts sound for all angles at such conditions. It does not correctly predict hot jets, jets of strong 
asymmetry, or jets with supersonic (relative to ambient speed of sound) exit velocities. 
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10. Appendix A. — Numerical Solution to Compressible Rayleigh Equation 

In a cylindrical coordinate system, the second order compressible Rayleigh operator governing 
function f m is written as (ref. 1 0) 


.. 1 p' M' 

f„ + 4 cos# f + 

m r p (l-Mcos#)_ 

co 2 2 w 2 co 1 2 /i 3cos 6 p' 1 

— (1-Mcos#) ^ y cos" 0 + — — — — (— M — M-M ) f m = 0, 

c r c x (l-Mcoso 1 ) p r 


where prime denotes a radial derivative. Having specified the jet profiles M = U(r)/c yi and 
c~ = ;/.H T(r), the numerical solution to (Al) starts with /„, = a m x m as r — >0 and continues to the jet boundary 
at r = R a where the function and its derivative are matched with the appropriate outgoing solution outside 
the jet. Constant a m thus becomes 


a , r k (J) - J'„ '(2-)] 

"ksmefJRJH*»(x)-dfJRJ/drH"'(x) 

{ 1, m = 0 

S m = ) 0 

[ 2 , m > 1 


(A2) 


In equation (A2), prime as appearing on the cylindrical Bessel and Hankel functions points to a 
derivative with respect to the argument j = kR a sin 6. 
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11. Appendix B. — Doppler Factor Singularity 

When a compact source is moving towards an observer at the speed of sound, equation (13) shows 
that the observer frequency as M c cos#->l. We may introduce a non-compact source in 

equation (6) by allowing the axial source position to be a complex number, with a small imaginary part 

Xj — ^ Xj (1 + i ccM c \ cc — ^0 (B 1) 

Mathematically, real parameter a allows a slight perturbation in the vicinity of the actual source 
location xj\ Equation (1 1) is written as 

_ 2n § {_ —[1 - M c {\ + i aM c ) cos 0\ + — } , a^O (B2) 

. u„ u „ 


Note that when the factor multiplying yq in the exponent of the above integrand is not zero, small 
parameter a should be selected appropriately to let the integrand converge at large y\ 


... . (O-CO 

• a -i(kcos0+—- — )y, 

e*'e Uc dy ] = 


/r-/r 


(P + -ix)(P -ix) la = a“-> 0 “,y 1 >0 


a = a + — > 0 + ,y 1 < 0 


and j3 = akM c cos#, j = Axos#-!- 


co -co 


Thus, the integral equation (B3) approaches zero as a—> 0, The delta function in 

equation (B2) is subsequently placed in equation (10). The integration with respect to co requires 

co s = 4(1 - M c cos 9)-iaM 2 c cos<9], as a — > 0. (B4) 

Since frequencies are real numbers, we may now define the modified Doppler factor as the magnitude 
of the above complex frequency 

co s = <»[(!- M c cos6*) 2 +(aM f ) 2 ] 1/2 , as M c cos#— >1 (B5) 


(1 - M c cos 0) — > [(1 - M c cos 0) 1 +a 2 M c 2 ] 1/2 (B6) 
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12. Appendix C. — Far-Field Spectral Density 

The acoustic pressure fluctuation at the observer point x and time t resulting from source Q emitting at point y 
and time t\ is 

P(x,t ) If ,t',y ,t\)Qiy ,ti)dt x dy , (Cl) 


where 


A +CO 

G(x,t;y,t x ) = — J G(x,y,co)e~ ,a>(t h) dco , 


(C2) 


-00 

and G(x,y, co) is an appropriate frequency-domain Green’s function. The sound spectral density is a Fourier 
transform of the auto-correlation function 


+oo +oo 

p 2 (x,co) = J e mT d t J p(x,t ) p(x,t + r)dt. 

-oo —oo 


(C3) 


Upon placing equation (C 1 ) in equation (C3) we fmd 

+oo +co 

p 2 (x, co) = J e U0T d r J dt\ J J G(x,t\ y, , t l )Q(y l , t l )dt x dy l 

-CO -CO [/j V, 


J J G(x,t + r;y 2 ,t 2 )Q(y 2 ,t 2 )dt 2 dy 2 

h h 


Equation (C2) is now placed in the preceding expression to obtain 


+oo +oo 


p 2 (x, co) = (— ) 2 J e im d r J dt^ J J J G(x, y 1 ,a) 1 )e m ' ( ' ,l) Q(v l , t x )da> l dt l dy 1 x 


-oo -co t x y { 


,y 2 ,o} 2 )e 


ieo 2 (t+T-t 2 


Q(y 2 ,t 2 )dco 2 dt 2 dy 2 k 


(C4) 


\co 2 t 2 y 2 


The integration with respect to t gives 2n:i%0)]+0h), which is followed by an integration with respect to ay, 
resulting in 


1 +00 

p 2 {x, co) = — J e“ or d x J J J J J G{x, y x ,co x )G{x,y 2 ,-co x )Q(y x , t x )Q(y 2 , t 2 )e m ' h 

-co ffl, V\ t 2 y 2 

e l( h) d x dt x dt 2 dy x dy 2 . 

The integration with respect to rgives a delta function 

r e ii(a+(0 ' )T dT = 27xS(co + co x ). 

J— 00 
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Subsequent integration with respect to oj\ results in 

p 2 (x,co) = ,y v -o))G(x,y 2 ,co)Q(y v t l )Q(y 2 ,t 2 )e , ® ( ' 1 h) dt l dt 2 dy l dy 2 . (C5) 

h Pi h y 2 

We now require that source function (f be a stationary function, i.e. Q = 0 when /| T , for some large time T. 

For a stationary source, the average value 

I Q(y, tf~- lim J I Q(y, tfdt (C6) 

21 T J T 

remains finite. Let 

(C7) 


1 +r 

q(y,co,T ) = — J Q(y,t)e uo ‘dt 


We define a time-domain cross-correlation 

1 +r 

i?(v,|,r) = 0‘(v 1 ,O0(v 2 ^ + z ') = lim — [£?XviT)0(v 2 T + r )^> ( C8 ) 

t^2T j t 

where c and r denote the spatial and temporal separations between points v, = y — 2; / 2, y 2 = y + g ! 2 and 
* indicates a complex conjugate. The Fourier transform of the cross-correlation function is given as 

+co 

Q n (y,},co) = J R(y,4, T)e UOT dr , (C9) 

-oo 

and is related to the spectral densities at points y, and y 2 

1 . 1 

lim — q (y l ,(Dj)q{y 2 ,6),T) = ——Q n (y,Z,G)). (CIO) 

r ->°°2 T (2ft) 

Expression (CIO) represents the spectral density for a two-point space-time correlation. Using the above 
equations in (C5) and noting that for a stationary process, the spectral intensity in now divided by 2 T we find 

— (2 jA 1 _ 

p 2 (x,co) = ^ J J G(x,y l -a>)G(x,y 2 , 0 })q\y x ,a>,T)q(y 2 ,a>,T)dy l dy 2 , (Cll) 

Pl?2 


or 

p 2 (x, co) = J J G* (x, y - 1 / 2 ,co)G(x, y +| / 2 ,co)Q n (y,g, co)d}dy. (C12) 

y 1 

Upon using equation (C9) in equation (C12), the spectral density per unit volume of turbulence at y is given as 
in equation (15). 


NASA/TM— 2005-213827 


22 



13. Appendix D. — Hase Variation of the Green’s Function 

The phase factor of a high-frequency Green’s function for a convecting source, in the notation of 
reference 5 is 


exp {ik[R + y((p- (p s ) + g n (r) - R sin 2 #]} , 
k 


(Dl) 


where R denotes a source to observer distance, which in the present notation should be replaced with 
(x-y ) , and 

r 

Sn( r ) = \QX r ) dr ’ ( D2 ) 

r s 


and r s is the zero crossing of the shielding function Q„{f), 


1-H/cos 

rk 


Q Z n (r) = (- 


c / c„ 


QM = o 


(D3) 


The phase variation of the GF between source points y x and y 2 becomes 

exp{/A'[A^j cos 0 + — A<p s +Ag n (/"')]} , 
k 


(D4) 


and 


'2 '2 

A ?„ (r s ) = J Qn ( r ) dr = ? sin e + J y/ n (r, G)dr 


(D5) 


The second term on the right hand side of (D5) signifies the effect of the refraction on phase. To the 
first order of approximation, this term may be neglected compared to the first term, which reflects the 
variations in phase due to source location only. As k— >co, equation (D4) becomes 

exp{ik[A<^ x cos 6* + A r s sin#]} , (D6) 


which is recognized as exp(7A . g) . 


NASA/TM— 2005-213827 


23 




14. Appendix E. — Space-Time Correlation Models 


When confronted with selecting the functional forms of the space-time correlation of velocity in a jet, 
several forms can be considered. Historically, these have been chosen partly based upon mathematical 
convenience, partly on data, and partly on compatibility with equations of motion. On the mathematical 
convenience, a model which has independent functions of space and time (in a moving frame) has often 
been chosen. And among functional forms for these independent functions, the Gaussian function has 
often been utilized for its ease of use in subsequent integrals. The Gaussian form also has the property of 
a zero derivative at the origin, making its Fourier transform integrable. 

Measurements made with hot-wire probes, shown here from reference 17, and supported by PIV 
measurements (ref. 16), show that a Gaussian function is not particularly well-suited for describing the 
form of the space-time correlation near the origin. These data were acquired in a Mach = 0.50 cold jet 
flow, and show the behavior of the correlation function in the axial and radial separations. Figures El and 
E2 show the correlation of the axial and radial components, respectively, of velocity at three different 
locations in the jet: x/Dj =4, r/Dj =0.50 (lipline); x/D, =8, r/Dj =0.50; and x/Dj =20, r/D, = 1.25. In each 
figure, the upper set of carpet plots show the correlation as a function of time delay and axial separation, 
while the lower plots show the correlation as a function of time delay and radial separation. There are 
interesting subtle variations in the correlations at the three locations in the jet, but they all show a very 
strong peaked behavior near the origin. Future versions may benefit from turbulence models which pick 
up the more subtle behaviors of these correlation functions, but at this time single-exponent models, such 
as equation 23, are better supported by data than Gaussian forms. 
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15. Appendix F. — Non-Compactness Factor 


When wavelength of the acoustic disturbances is small relative to the correlation length scale 
(i.e., A < £), it is shown (ref. 1) that the axial correlation coefficient /mi (defined in eq. 25) consists of a 
non-compactness factor/V(AT ) at high frequency. For example, using an exponential spatial function 
(see eq. 23) in the two-point correlation function of Batchelor, the axial correlation component R n 
becomes 


R,m=ui[(f+\j.ns, 


(FI) 


m = e-« n 

(F2) 


Note that for brevity, £is used here as the separation vector in a moving frame. Using (F2) in (FI) 
results in 


R n {%) = ul\\ 


n 

2QL 




-n£lt 


(F3) 


Upon placing (D3) into equation (25) we find 

I nn {y,co s ) = 2 H{co s )\ (F4) 

1 

The above integration my be performed in a spherical coordinate system 
£ = (^cosa, ^ sin a cos sin a sin ^)) and by selecting dummy variable d, such that d, x aligns with 
direction of radiation k 

co 2k 71 

I un = 2(u x ) 2 £ 3 H(g> s ) J J J (1-— £sin 2 a) 2 exp(-2^-/(A:^cosa)(^ 2 sinadadcfxi^) (F5) 

£=0 <f= 0 a=0 ^ 


This leads to equation (27) with the non-compactness factor given as 

U 


jr 

N(k£)= 20(-) 5 [3tan , (— )-2- 


5(— ) 2 + 12 


n 


kt 


'2k 


K 


(— ) 2 +4 

V 7 t 


], for /(£) = e 


-n£ll 


(F6) 


The non-compactness factor thus becomes a function of the spatial form of the correlation. For 
example, when a Gaussian function f(^) = exp(-;r<^ 2 / 1 2 ) is placed in (FI), equation (27) should be 
calculated from 


N(kl) = exp( 


2/>2 


k z £ 


8 n 


for 


m = e 


Tcfn 2 


(F7) 
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17. Figures 



Figure 1 . — Spectral shape function at 90° with 
temporal function defined as: exp(-|r/rj), 
solid line; exp{-(0.4 2 + r 2 /r 2 ) 0 50 }, dash-dot; 
exp (-r 2 /r 2 ), dashed line. 



Figure 2(a). — Directivity D(x,x\ to) due to a stationary ring 
source in a Mach 0.9 unheated jet at St = 0.10, 
x‘ !Dj= 1.0 and r/Dy 0.0 (solid line); 0.50 (dashed line); 
1.0 (dash-dot). 
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Figure 3(a). — Directivity D(x,x',co) due to a 
stationary ring source in a Mach 0.9 unheated 
jet at St = 0.10, x[ /Dj= 7.5 and r/Df. 0.0 (solid 
line); 0.50 (dashed line); 2.0 (dash-dot). 



Figure 3(b). — As figure 3(a) but for St = 0.25. 



Figure 3(c). — As figure 3(a) but for St = 1.0. 
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Figure 3(d). — As figure 3(a) but for St = 5.0. 



Figure 4. — Directivity D(x,x s ,a>) due to a convecting 
ring source in a Mach 0.9 unheated jet at observer 
St = 0.25, x'l /Dj= 1.0 and r/D/. 0.0 (solid line); 
0.50 (dashed line); 1.0 (dash-dot). 



Figure 5. — As figure 4 but for a ring source at 
x\IDj=1.5 and A/D y. 0.0 (solid line); 

0.50 (dashed line); 2.0 (dash-dot). 
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. — Spectra in Mach 0.51, T r =1.0 (SP03) jet at indicated inlet angles and at R/Dj= 100. 
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Figure 7(a). — Spectra in Mach 0.98, T r =1.0 (SP07) jet at indicated inlet angles and at R/Dj= 100. 
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Figure 7(b). — Comparison between calculated spectra and data at indicated inlet angles at 
r/Dj= 100. (a) AO =1.1 85, (b) AO =1.40, (c)AO=1.50, (d)AO= 1.66,(e)AO= 1-80. 



Figure 8(a). — Directivity A (x,x\a>) due to a 
stationary ring source in Mach 0.9, T r = 2.78 jet 
at St = 0.10, x[IDj=1.5 and r/Df. 0.0 
(solid line); 0.50 (dashed line); 2.0 (dash-dot). 
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Figure 12(a). — Distribution of turbulent kinetic energy (SP03). 

Top half, RANS prediction; bottom half, PIV measurements. 
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igure 12(b). — Distribution of turbulent kinetic energy (SP07). 

Top half, RANS prediction; bottom half, PIV measurements. 
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Figure 12(c). — Distribution of turbulent kinetic energy (SP49). 
Top half, RANS prediction; bottom half, PIV measurements. 
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Figure El. — Space-time correlation of u x u x . 



Figure E2. — Space-time correlation of u 2 u 2 ' . 


A comparison of the measurements with model forms using Gaussian 
and exponential functions is shown in reference 28. 
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